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We present evidence for D^-D^^ mixing in D" — > K'^-k^ decays from 384 fb~^ of e^e~ coUiding- 
beam data recorded near ^/s = 10.6 GeV with the BABMI detector at the PEP-II storage rings 
at SLAG. We find the mixing parameters x'^ = [-0.22 ±0.30 (stat.) ±0.21 (syst.)] x 10"^ and 
y' = [9.7±4.4 (stat.)±3.1 (syst.)] x 10"^, and a correlation between them of —0.94. This result 
is inconsistent with the no-mixing hypothesis with a significance of 3.9 standard deviations. 
We measure Ru, the ratio of doubly Cabibbo-suppressed to Cabibbo-favored decay rates, to 
be [0.303 ± 0.016 (stat.) ± 0.010 (syst.)]%. We find no evidence for CP violation. 

1 Introduction 

The and-D*^ mesons are flavor eigenstates which are invariant in strong interactions, but are 
subject to electroweak interactions that permit an initial flavor eigenstate to evolve into a time- 
dependent mixture of and D^. In the Standard Model (SM), such oscillations proceed through 
both short-distance and non-perturbative long-distance amplitudes. The expected mixing rate 
mediated by down-type quark box diagrams and di-penguin^ diagrams is O{10~^ — lO"^*^), 
while the predicted range for non-perturbative long-distance contributions is approximately 
bounded by the box diagram rate and the current experimental sensitivity of 0(10^'*). New 
physics predictions span the same large range, We present evidence for L) mixing consistent 
with these expectations and with previous experimental limits.' ^ 6 7 8 9,10,1 1 ^^^^ compare 
and samples separately, and find no evidence for CP violation. 

The mixing rate is characterized using the right-sign (RS), Cabibbo-favored (CF) decay 
K~Tr~^ and the wrong-sign (WS) decay K~^tt~ . The WS final state can be produced 

either through a doubly Cabibbo-suppressed (DCS) tree-level decay or through mixing followed 
by a CF decay. The DCS decay has a small rate Rd of order tan^ 9c ~ 0.3% relative to CF 
decay, where 9c is the Cabibbo angle. We distinguish and by their production in the 
decay D*~^ t:'^D^ . In RS decays, the 7r+ and kaon have oppositely signed charges, while in 
WS decays the charge signs are the same. The time dependence of the WS decay rate is used 
to separate DCS from mixed decays. 

Charm mixing is generally characterized by two dimensionless parameters, x = Am/F and 
y = AF/2F, where Am = m2 — mi {AT = F2 — Fi) is the mass (width) difference between the 
two neutral D mass eigenstates and F is the average width. If either x or y is non-zero, then D 
mixing will occur. We approximate the time dependence of the WS decay of a meson produced 



"The use of charge-conjugate modes is implied unless otherwise noted. 



as a at time t = in the limit of small mixing (|x|, \y\ <^ 1) and CP conservation as^^^l 

^ocRo + y^y' Tt + '^^^±y^{Ttf , (1) 

where x' = x cos Sx-K+y sin 5k-w, y' = —x sin 5Kn+y cos 6Kn, and 6ktt is the strong phase between 
the DCS and CF amplitudes. We study both CP-conserving and CP-violating cases. For the 
CP-conserving case, we fit for the parameters Rd, x''^, and y'. To search for CP violation, we 
apply Eq. [T] to the and samples separately, fitting for the parameters , x'"^^ , y'^ for 
P»°(+) and DO(-) decays. 

2 Event Selection and Analysis 

We use 384 fb"^ of e+e- coUiding-b eam data recorded near -v/s = 10.6 GeV with the Babar 
detector^l^ at the PEP-II asymmetric-energy storage rings. We initially select signal candidates 
by combining oppositely-charged tracks identified as -fC or vr using a likelihood-based particle 
identification algorithm, requiring the K^ir^ invariant mass 1.81 < ruKn < 1-92 GeV/c^ and 
e~^e~ center-of-mass frame (CM) momentum p^Q > 2.5 GeV/c. We require the vr^ to have 
laboratory momentum p > 0.1 GeV/c and CM momentum p* < 0.45 GeV/c. 

To obtain the proper decay time t and its error at for each candidate, we refit the 
daughter tracks and the vr^, constraining the daughters to originate from a common vertex 
while simultaneously requiring the and 7r+ to originate from a common vertex constrained 
by the position and size of the e~^e~ interaction region. We require a refit probability 
Pix^) > 0.1%, M^Kmrf) — rriKn mass difference 0.14 < Am < 0.16 GeV/c^, proper decay time 
—2 < t < 4 PS and proper decay time error at < 0.5 ps. The nominal mean proper lifetime is 
~ 0.410 ps^^and the most probable value of at for signal events is 0.16 ps. If there are multiple 
signal candidates with overlapping tracks within an event, we retain only the candidate with the 
highest P(x^). After applying all selection criteria, we retain approximately 1,229,000 RS and 
64,000 WS signal candidates. To avoid potential bias, the complete event selection and analysis 
procedures were finalized prior to examining the mixing results. 

The mixing parameters are determined using an unbinned, extended maximum-likelihood 
fit to the RS and WS data samples over the four observables rriK-K, Am, t, and at- The fit 
is performed in several stages. First, the shapes of the RS and WS signal and background 
probability density functions (PDFs) are determined from an initial 2-d fit to {mx7r,Am}. 
These shapes are then fixed in subsequent fits. Next, the proper-time resolution function 
and lifetime are determined from a fit to the RS data using {rriK-K, Am} to separate the signal 
and background components. Finally, the WS data is analyzed using three different fit models. 
The first model assumes both CP conservation and the absence of mixing, the second model 
allows mixing but no CP violation, and the third model allows mixing and CP violation. 

The RS and WS {mx-K-, Am} distributions are described by four components: signal, random 
vr^, misreconstructed and combinatorial background. The signal component has a character- 
istic peak in both mx-w and Am. The random vr^ component models reconstructed decays 
combined with a random slow pion and has the same shape in mx-w as signal events, but does 
not peak in Am. Misreconstructed events have one or more of the decay products either 
not reconstructed or reconstructed with the wrong particle hypothesis. They peak in Am, but 
not in mxn- For RS events, most of these are semileptonic decays. For WS events, the main 
contribution is RS K~tt^ decays where the K~ and the tt"*" are misidentified as tt~ and 

K^, respectively. Combinatorial background events comprise the remainder of events and do 
not exhibit any peaking structure in mx-K or Am. 

We fit the RS and WS data samples simultaneously to determine the PDF parameters 
describing the signal and random vr^ event class shapes shared between RS and WS datasets. 



We find 1, 141, 500 ± 1, 200 RS signal events and 4, 030 ± 90 WS signal events. The dominant 
background component is the random vr^ background. Projections of the WS data and fit are 
shown in Fig. [TJ 
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Figure 1: a) niK-K for WS candidates with 0.1445 < Am < 0.1465 GeV/c^, and b) Am for WS candidates with 
1.843 < niK^ < 1.883 GeV/c^. The fitted PDFs are overlaid. 

The measured proper-time distribution for the RS signal is described by an exponential 
function convolved with a resolution function whose parameters are determined by the fit to the 
data. The resolution function is the sum of three Gaussians with widths proportional to the 
estimated event-by-event proper-time uncertainty at- The random vr^ background is described 
by the same proper-time distribution as signal events, since the slow pion has little weight in the 
vertex fit. The proper-time distribution of the combinatorial background is described by a sum 
of two Gaussians, one of which has a power-law tail to account for a small long-lived component. 
The combinatorial background and real decays have different at distributions, as determined 
from data using a background-subtraction technique^^^ based on the {niKn, ^n^} fit- 

The fit to the RS proper-time distribution is performed over all events in the full mKir and 
Am region. The PDFs for signal and background in rnxn and Am are used in the proper-time 
fit with all parameters fixed to their previously determined values. The fitted lifetime is 
found to be consistent with the world-average lifetime. 

The measured proper-time distribution for the WS signal is modeled by Eq. [T] convolved 
with the resolution function determined in the RS proper-time fit. The random vr^^ and mis- 
reconstructed backgrounds are described by the RS signal proper-time distribution since 
they are real decays. The proper-time distribution for WS data is shown in Fig. [2l The fit 
results with and without mixing are shown as the overlaid curves. The fit with mixing provides 
a substantially better description of the data than the fit with no mixing. The significance of 
the mixing signal is evaluated based on the change in negative log likelihood with respect to 
the minimum. Figure [3] shows confidence-level (CL) contours calculated from the change in 
log likelihood (— 2Aln£) in two dimensions {x'"^ and y') with systematic uncertainties included. 
The likelihood maximum is at the unphysical value of = -2.2 X 10"^ and y' = 9.7 x lO'^. 
The value of — 2Aln£ at the most likely point in the physically allowed region (x'^ = and 
y' = 6.4 X 10~^) is 0.7 units. The value of — 2Aln£ for no-mixing is 23.9 units. Including the 
systematic uncertainties, this corresponds to a significance equivalent to 3.9 standard deviations 
(1 — CL = 1 X 10~^) and thus constitutes evidence for mixing. The fitted values of the mixing 



parameters and Rd, along with errors, are listed in Table [H The correlation coefficient between 
the x'^ and y' parameters is —0.94. 




Figure 2: a) Projections of the proper-time distribution of combined D° and WS candidates and fit result 
integrated over tiie signal region 1.843 < tuktv < 1-883 GeV/c^ and 0.1445 < Am < 0.1465 GeV/c^. The result 
of the fit allowing (not allowing) mixing but not CP violation is overlaid as a solid (dashed) curve, b) The points 
represent the difference between the data and the no- mixing fit. The solid curve shows the difference between fits 

with and without mixing. 




Figure 3: The central value (point) and confidence-level (CL) contours for 1 — CL — 0.317 (Icr), 4.55 x 10 ^ (2o"), 
2.70 X 10"^ (3cr), 6.33 x 10"^ (4cr) and 5.73 x 10"^ {5a), calculated from the change in the value of -21n£ 
compared with its value at the minimum. Systematic uncertainties are included. The no-mixing point is shown 

as a plus sign (-I-) . 



Allowing for the possibility of CP violation, we calculate the values of Rd = \J Rj^Rj^ and 
Ad = {Rd - Rd)/{Rd + ^d) listed in Table H from the fitted values. The best fit points 
(x'^^,?/'^) shown in Tabled] are more than three standard deviations away from the no-mixing 
hypothesis. The shapes of the (x'^^,y'^) CL contours are similar to those shown in Fig. [3l All 
cross-checks indicate that the close agreement between the separate and fit results is 
coincidental. 

As a cross-check of the mixing signal, we perform independent {mx-K^ ^fn{ fits with no 
shared parameters for intervals in proper time selected to have approximately equal numbers 



Table 1: Results from the different fits. The first uncertainty listed is statistical and the second systematic. 



Fit type 


Parameter 


Fit Results (/lO"^) 


No CP violation or mixing 


Rd 


3.53 ± 0.08 ± 0.04 


No CP violation 


Rd 


3.03 ±0.16 ±0.10 




X 


-0.22 ±0.30 ±0.21 




y' 


9.7 ±4.4 ±3.1 


CP violation allowed 


Rd 


3.03 ±0.16 ±0.10 




Ad 


-21 ± 52 ± 15 




x' ^ 


-0.24 ± 0.43 ± 0.30 






9.8 ±6.4 ±4.5 




x'^~ 


-0.20 ±0.41 ±0.29 




y'- 


9.6 ±6.1 ±4.3 



of RS candidates. Figure H] shows the resulting fitted WS branching fractions growing with 
increasing proper time. The slope of a linear fit to the data points is consistent with the 
measured mixing parameters and inconsistent with the no-mixing hypothesis. 

0.45 1 



0.4 



0.35- 



0.3 



-2-10 1 2 

t(ps) 

Figure 4: The WS branching fractions from independent niK-n-, Am fits to slices in measured proper time (points). 
The dashed line shows the expected wrong-sign rate as determined from the mixing fit shown in Fig. [5] The 
with respect to expectation from the mixing fit is 1.5; for the no-mixing hypothesis (a constant WS rate), the 

is 24.0. 



We validated the fitting procedure on simulated data samples using both MC samples with 
the full detector simulation and large parameterized MC samples. In all cases we found the fit to 
be unbiased. As a further cross-check, we performed a fit to the RS data proper-time distribution 
allowing for mixing in the signal component; the fitted values of the mixing parameters are 
consistent with no mixing. In addition we found the staged fitting approach to give the same 
solution and confidence regions as a simultaneous fit in which all parameters are allowed to vary. 

In evaluating systematic uncertainties in i?/) and the mixing parameters we considered vari- 
ations in the fit model and in the selection criteria. We also considered alternative forms of 
the raK-K^ Sm, proper time, and at PDFs. We varied the t and at requirements. In ad- 
dition, we considered variations that keep or reject all D*^ candidates sharing tracks with 
other candidates. For each source of systematic error, we compute the significance sf = 



lnC{x 



y') 



/2.3, where {x' , y') are the parameters obtained from the stan- 

the parameters from the fit including the i^^ systematic variation, and C the 

i confidence level for 2 degrees of free- 



dard fit, (x'^,y-) 

likelihood of the standard fit. The factor 2.3 is the 



dom. To estimate the significance of our results in {x' , y'), we reduce — 2Aln>C by a factor of 
1 + Esf = 1.3 to account for systematic errors. The largest contribution to this factor, 0.06, 
is due to uncertainty in modeling the long decay time component from other D decays in the 
signal region. The second largest component, 0.05, is due to the presence of a non-zero mean 
in the proper time signal resolution PDF. The mean value is determined in the RS proper time 
fit to be 3.6 fs and is due to small misalignments in the detector. The error of 15 x 10^^ on Ajy 
is primarily due to uncertainties in modeling the differences between and K~ absorption in 
the detector. 

We have presented evidence for D^-D^ mixing. Our result is inconsistent with the no-mixing 

hypothesis at a significance of 3.9 standard deviations. We measure y' = [9.7 it 4.4 (stat.) it 
3.1 (syst.)] X 10^^, while x'^ is consistent with zero. We find no evidence for CP violation and 
measure Rd to be [0.303 ± 0.016 (stat.) ± 0.010 (syst.)]%. The result is consistent with SM 
estimates for mixing. 
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